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Abstract 

The halo dark matter (DM) can be captured by the Sun if its final velocity after the collision 
with a nucleus in the Sun is less than the escape velocity. For self-interacting DM (SIDM), we 
show that the number of DM trapped inside the Sun still provides significant signals even if the 
DM-nucleon cross section is negligible. We consider a SIDM model where U{1) gauge symmetry is 
introduced to account for the DM self-interaction. Such a model naturally leads to isospin violating 
DM-nucleon interaction, although isospin symmetric interaction is still allowed as a special case. 
We show that the detection of neutrino signature from DM annihilation in the Sun can probe those 
SIDM parameter ranges not reachable by direct the detection. Those parameter ranges are either 
the region with a very small m x or the region opened up due to isospin violations. 
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I. INTRODUCTION 


Single component and collisionless cold dark matter (CCDM), which is treated as the 
standard dark matter (DM) candidate in ACDM model, accounts for the cosmological data 
from Cosmic Microwave Background (CMB), Big Bang Nucleosynthesis (BBN), and the 
large scale structure. However, the above DM properties lead to some controversies between 
the N-body simulation and astrophysical observations on small scale structures with non¬ 
linear DM-dominated systems. A cusp structure is formed by gravitationally attracting a 
large number of DM in the center regions of galaxies and the zero dissipative dynamics [Ij. 
On the other hand, observations indicate that DM distributes in a much more flat profile in 
the center regions m- This is called the cusp/core problem [4j. There exists other problem 
with CCDM concerning the size of subhalos, which are the hosts of the satellite galaxies 
surrounding the Milky Way (MW) halo. The dispersion velocities from these galaxy rotation 
curves reflect the size of their host subhalos. It is observed that a discrepancy in the subhalo 
size exists between the CCDM-only simulation and the observations [SI IB] - About 0(10) 
most massive subhalos generated from the N-body simulation are too massive in the MW 
halo (with circular velocity larger than 30 km/s) whereas the maximum circular velocities 
of MW dwarf spheroidals are less than 25 km/s. This is referred to the “Too-big-to-fail” 
problem indicating no satellite galaxy hosted by such massive subhalo is found. 

It is possible that these puzzles may be due to insufficient knowledge of the baryonic pro¬ 
cesses such as the supernova feedback and photonization [MS]- On the other hand, these 
conflicts may also hint nontrivial features of DM such as self-interacting DM (SIDM) [IT] . 
It was then noted that the SIDM with a constant cross section cannot account for observed 
cllipticities in clusters m H6 j and the survivability of subhalos H3- In the very short 
mean free path limit, it produces even more cuspy profiles [T8li2T] . However, some recent 
simulations have shown that the SIDM velocity-independent cross sections in the range 
0.1 cm 2 /g < o' 7 ’/m x < 10 cm 2 /g ( m x denoting the DM mass) can resolve the cusp/core and 
Too-big-to-fail problems on dwarf scales [22l [2o] . although it was pointed out in Ref. [23] 
that orlvri x ~ 0.1 cm 2 /g is too small to account for the population of massive subhalos, and 
the study using bullet cluster IE 0657-56 constrains to be less than 1.25 cm 2 /g [26]. 

Furthermore, other investigations of SIDM with characteristic velocity-dependent cross sec¬ 
tions provide a broader cross section range, 0.1 cm 2 /g < o^jm x < 50 cm 2 /g, for alleviating 
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the above-mentioned puzzles [27H3T] 1 . From particle physics point of view, DM exchanging 
light mediators would generate an attractive self-interaction, which enhances the annihila¬ 
tion cross section by Sommerheld effect and resonance scattering (33] EE]- The mediator 0 
often plays as the messenger between the visible matter and dark matter. This mediator can 
be scalar, pseudoscalar, vector, or axial-vector particles [36]. The scalar and vector interac¬ 
tions with nucleon only generate spin-independent (SI) cross section (in the non-relativistic 
limit), while spin-dependent cross section can be generated via exchanging pseudoscalar 
or axial-vector particles. The relatively large SIDM cross section requires a MeV scale 0 
satisfying m 


(j ~ 

u xx 


7.6 x 10 -24 cm 2 


a, 


. m 4> y*( m x y 

0.01/ V 30 MeV / vGeV/ 


( 1 ) 


where a x is the dark fine structure constant representing the coupling strength between DM 
X and the mediator 0. With m $ around the MeV scale, it is possible for 0 to decay into light 
standard model particles. In such a case, the lifetime of 0 should be less than 1 sec to satisfy 
the BBN constraint. 2 Therefore, the light mediator 0 can produce both direct and indirect 
DM signatures. The former signature is generated when DM scatters with the nuclei by 
exchanging 0. The latter signature is generated by the annihilation of DMs into a pair of 0, 
which subsequently decay into SM particles. The mediator 0 can be a hidden 17(1) gauge 
boson or scalar, which mixes with photon, Z boson or Higgs boson to bridge between DM and 
the visible sectors. It should be noted that Eq. (JTJ) is valid only at the perturbative regime 
(a x m x /< 1), while the effects of dark force on the halo structure for the full range 
of parameters were investigated in Ref. jdOj. The characteristic velocity-dependent cross 
section at different scales for resolving all the conflicts can be realized. Furthermore, DM 
direct detection experiments for probing such SIDM scenarios were proposed in Ref. [3Ej (see 
also, for example, Ref. [39] for the asymmetric DM scenario). The predicted direct detection 
cross section is within the reach of next-generation experiments such as XENON1T HI and 
SuperCDMS [31 j. 

In this paper, we shall consider symmetric fermionic DM and show that the annihilation 
signatures from the trapped DMs inside the Sun can provide a complementary test to SIDM 


1 The thorough study on cosmology and structure formation in the context of mirror photon is investigated 
in Ref. [33]. Some early attempts to study the astrophysical effects of mirror matters can also be found, 
for example, in Ref. [33] , 

2 If DM is kinematically decoupled from the dark radiation at late time, it would leave signals in CMB EZ]. 
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scenarios. It is well known that the DM-nucleon scattering cross section relevant to the 
indirect DM signature from the Sun is identical to the cross section relevant to DM direct 
detection experiments. For the former case, such a cross section leads to the DM capture 
if the velocity of the final-state DM is less than the escape velocity of the Sun. As the 
number of captured DMs increases to a significant level, the rate of DM annihilations into 
SM particles could become detectable H2H5H. The total number of the captured DMs is 
determined by the above DM-nucleon scattering cross section and the DM mass. However, 
it has been shown in a model-independent approach that the DM accumulation time can 
be shortened in the Sun and the total number of trapped DM can also be increased if 
SIDM is considered [491I5T] . In Ref. [SO], the constraint ot < 0( 1CT 22 ) cm 2 is obtained 
for the xx —> W + W~/t + t~ annihilation mode using IceCube 79 string data [52] which 
probes the mass range m x > 20 GeV. The low-mass DM region to be probed by IceCube- 
PINGU is considered in Ref. |5T]. In this work we shall present a realization of the above 
model-independent results with a concrete SIDM model. We work out the SIDM model 
prediction on the DM-induced neutrino flux from the Sun and the sensitivities of future 
neutrino telescopes to such a signal and the related model parameters. We compare these 
sensitivities to constraints set by current direct detection experiments. We will show that 
the searches for direct and indirect DM signals give complementary tests to the parameter 
space of SIDM model. The complementarity of direct and indirect searches arises from two 
reasons. First, the indirect search can look for very light DM which is not yet probed by 
direct detections. Second, the direct detection sensitivity could be significantly worsen by 
isospin violation in DM-nucleon couplings. On the other hand, the DM-induced neutrino 
flux from the Sun can be insensitive to DM-nucleus scattering cross section so long as DM 
self-interaction is strong enough j49j, [51]. 

Our paper is organized as follows. In Section II we briefly introduce the SIDM models 
where the light mediator is either a gauge boson or a scalar particle. The constraints on 
these models are summarized. Section III shows that the DM accumulation inside the Sun 
can be significantly enhanced by DM self-interaction. In Section IV, we demonstrate that 
IceCube/PINGU experiment can probe SIDM parameters in a way complementary to the 
direct detection experiment. We present our conclusions in Section V. 
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II. SCATTERING VIA LIGHT MEDIATOR EXCHANGE 


Two simplest ways to construct the SIDM model with the feature of Sommcrfield en¬ 
hancement are to assume a light vector boson or a light scalar boson as the dark force 
mediator 3 . For the vector mediator case, DM carries a charge en under a hidden Abelian 
U( l) x gauge symmetry and 0^ is the corresponding gauge boson. The hidden U( l) x gauge 
boson has no ordinary SM interaction except the one induced by kinetic mixing with the SM 
photon [53] ■ Since we expect the mass of dark mediator 0 to be around MeV scale, there 
must exist certain mechanism to generate a massive U(l) x gauge boson. It can be certain 
spontaneous symmetry breaking pattern in the dark sector or based upon the Stueckelberg 
mechanism [HU 55] to ensure the UV completion of the theory. Here we take the phenomeno¬ 
logical approach and neglect the details of model construction. The 0y and (j)Z mixings lead 
to the following couplings between 0^ and SM currents: 


£mixing/vector ( T 6z 

Cw 




with 


JS, = Ee//y/ 4 c = E-tv 


'3/ 


1-75 


QfS 


w 


f 


( 2 ) 


( 3 ) 


/ / 
representing electromagnetic current and weak neutral current respectively. Here / is the 
SM fermion and Qf refers to its electric charge, g 2 is the SU(2) L gauge coupling, and 
cw{sw ) = cos6 l vi/(sin6 l w). The parameters e 7 and are originated from kinetic and mass 
mixings between 0 M and gauge bosons such as 


l-'mixing, [/(l) T ' 


With DM x a Dirac fermion, the dark force is given by 


C 


DM, 


^vector 


= e D xYx$n- 


( 4 ) 


( 5 ) 


For the case of scalar 0, it is possible that 0 can mix with Higgs boson. The simplest 
scenario is that 0 is a real singlet scalar. The relevant terms of the Lagrangian are then 
given by 


£ 


DM, 


^scalar 


3 g,XX<P + a<t>\H\ 2 + b<p 2 \H\ 2 , 


( 6 ) 


In this paper, we analyze the case of velocity-independent cross sections. 
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where g s , a , and b are coupling constants and H is the SM Higgs doublet. As the Higgs 
boson develops the vacuum expectation value v ~ 246 GeV, the mass matrix between 0 
and H can mix via the a term. We can define the mixing parameter ~ av/m 2 h , thus the 
effective Lagrangian for the coupling of scalar 0 to the SM fermions is given by 


£ 


mixing/scalar 


E-^/ 7 . 

v 

f 


(7) 


A. Scattering with vector mediator 


The Feynman diagrams for scattering among two Dirac fermionic (anti-) DMs via the dark 
mediator 0 are depicted in Fig. [lj Due to Fermi statistics, the u and t channel amplitudes 
of XX scattering carry a relative minus sign such that two amplitudes cancel each other 
in the non-relativistic limit. The same cancellation occurs in XX scattering. Hence only 
XX scattering cross section is nonvanishing in the non-relativistic limit. Furthermore, the s 
channel amplitude of XX scattering is suppressed so long as m x . Hence we consider 

only t channel amplitude for XX scattering. The perturbative calculation gives the cross 
section 

9 

Tfl 

^xx~47ro7-f, (8) 

"V 

where a x = e|,/(47r) is the fine structure constant in the hidden 17(1) sector. Since only a xx 
is not suppressed for Dirac DM, the parameterization in Eq. (|TJ) should be understood as the 
one for a xx . In fact, for simplifying the notations, we henceforth denote the XX scattering 
cross section as a xx . 

The DM-nucleon scattering is shown in Fig.[2]where N stands for either proton or neutron. 
The parameter e^v is the strength of 0^-nucleon coupling in the unit of electric charge e. 
The SI cross section between DM and any nucleus with mass number A and proton number 
Z is then given by 


°xA 


167ia x a e 

ml 


\e v Z + e n (A — Z)] n xA — 


167ra x a e 

ml 


' £ p[ z + ri(A - Z)Y^ xA , (9) 


with a em the fine structure constant, [i x a = m x mA/ (m x + tua) the reduced mass for DM- 
nucleus system and 7] = £ n /£ P is the isospin violation parameter. The cross section for x 
is identical such that o x a = 0 x a- The effective 0^-nucleon couplings, e p and e n , can be 
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FIG. 1: The Feynman diagrams contributing to DM scatterings. The upper panel represents the t 
and u channel diagrams of XX scattering. The lower one represents the s and t channel diagrams 
of XX scattering. The XX cross section vanishes in the non-relativistic limit while XX scattering 
is dominated by the t channel diagram. The behavior of XX scattering is the same as that of XX 
scattering. See the main text for details. 
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FIG. 2: The Feynman diagram contributing to ,\-iV scattering. One simply replaces x by X ' n the 
diagram for x-IV scattering. 


expressed in terms of £ 7 and such that [38] 




£ £ 

£p — X + (1 4s 

4 S\yC\y 

2 \ ^ 
w) ~ 

£ 7 T 0.05 Ez , 

(10a) 

£n- . « 0.6 

4 S\yCw 



flOb) 

Thus, we have SI cross sections 




cr®p ra 1.5 x 10 -24 cm 2 e 2 ( 

' a x N 

\ ( m <t> V 4 

(11) 

,0.01. 

1 V30 MeV / 

for DM-proton scattering and 




a s ^ n a 5 x 1CT 25 cm 2 e\ ^ 

a x \ 

/ \ ~ 4 

(12) 

0.01/ 

V 30 MeV/ 


for DM-neutron scattering. 
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B. Scattering with scalar mediator 


The Feynman diagrams for DM-DM scattering with scalar exchange are identical to those 
depicted in Fig. [lj Tims we only replace a x by a s = g 2 /47r, and the mixing parameter e h 
via the Higgs boson are the same for both proton and neutron. Roughly one has 

£p,n ~ 3 x HT 3 ^. (13) 


The general DM-nucleon cross section thus obeys the isospin symmetry with 


a xN 


2 x 10 


-2y cm 2 e 2 , 


a s \ / \ 

h l(hOW V30MeV/ 


-4 


(14) 


However, to have (p decay before BBN, one requires Eh > 10 -5 . With such a lower bound 
for Eh , the resulting <j xp is so large that it is excluded by the current direct search for 
m x > 10 GeV [38]. In this work, we shall not continue discussing the scalar case, although 
some other scalar-exchange models might still be viable. It is important to note that the 
vector mediator case already contains many features of SIDM models. In addition, we 
are interested in searching for neutrino signals. Such type of signals disfavors the scalar 
mediator since the decay branching ratios of <p to neutrinos are generally negligible, unless 
?7i0 < 1 MeV. We note that the spin-dependent DM-nucleon cross section is suppressed by 
0(l/m x ) for the case of vector mediator. On the other hand, such suppression does not 
occur if <p is an axial-vector particle. In this paper, we shall only focus on spin-independent 
cross section. Hence we simplify cr®* as a xp from now on. 


III. DM ACCUMULATION IN THE SUN 


A. The evolution equation 


Since we consider the scenario that x and X are equally populated, the halo DM number 
density near the solar system can then be written as po = p x + p x = 0.3 GeV cnr 3 with 
Px = Px = 0-15 GeV cnr 3 . The DM in the Sun consists of two species, the DM and anti-DM 
with the numbers N x and N x respectively. The time evolutions of N x and N x are given by 


dN x 

dt 

dN x 

dt 


Cc - C e N x + C S N X - (C a + C se )N x N x , 
C c - C e N x + C,N X - (C a + C se )N x N x , 
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(15a) 

(15b) 






with C c the capture rate, C e the evaporation rate, C s the capture rate due to self-interaction, 
C se the self-interaction induced evaporation rate, and C a the annihilation rate. They are 
taken to be time-independent and have been fully discussed in Refs. (UHlHlEIj and references 
therein. 

The capture rates C c and C s depend on DM-nucleus scattering cross section g x a and 
DM-DM self-interaction cross section a xx , respectively. For a MeV range m^, it has been 
pointed out that [38, Ml g x a is sensitive to the momentum transfer q flowing into the (j) 
propagator shown in Fig. [2} We thus have 


a x A(q 2 ) 


m * -0 

K + 9 2 ) 2 


(16) 


where the magnitude of momentum transfer is given by q 2 = 2 m A E R with Er the nucleus 
recoil energy, and a^ A = a xA (q 2 = 0) is the cross section with zero momentum transfer. It 
is clearly seen that <j x a(q 2 ) is suppressed compared to cr° A . For estimating this suppression, 
we take E R = 2 n x A v x /mA with /j, x a the reduced mass for the DM-nucleus system. This 
Er corresponds to an 180° DM recoil angle after the collision in the center of momentum 
frame. We further take v x ~ 270 km/s, which is the DM velocity dispersion in the halo. 
Apparently, q 2 depends on both in a and m x . Thus, the capture rate with the momentum 
transfer suppression can be expressed as 


C, oc 


Px 


0.15 GeV/cm 3 y \ m 


( GeV\ / 270 km/s 


F Ai.m x ,i])a 0 X A- 


rnz 


mi + q A 


2\2 ' 


(17) 


where FA(m x , i]) is the product of various factors relevant to the chemical element A in 
the Sun, including the mass fraction, chemical element distribution, kinematic suppression, 
form factor, reduced mass and isospin violation effect. In Fig. [3j we compare the effect of 


momentum transfer suppression on C c for different m ?J with cr^ p fixed at 10 45 cm 2 . One can 
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see that the momentum transfer suppression is more severe for larger m x since q A 2 increases 
The above momentum transfer suppression also occurs in C s . The suppression factor can 


with m x . 


be taken from Eqs. (16) and (17). Thus 


C s (q 2 ) = C s (q 2 = 0). 


m. 


+ 9 


2^2 


(18) 


For DM self-capture, we have q = m x {vi — Vf) where Uj is the DM velocity before the 
scattering with the magnitude of the velocity given by u* = Jv^+v esc (r) 2 , while Vf is 
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T] = 1 , &xp = 10 ^ cm^ rj = - 0 . 7 , (r^ p = 10 ^ cm^ 




FIG. 3: The capture rates with (solids) and without momentum transfer (dashed). The cross 
section a® p = 10 -4,5 cm 2 for all results. Left and right are corresponding to ?y = 1 and —0.7. 


the DM velocity after the scattering with Vf = v esc (r )/\/2 on average ra- We take the 
approximation that Vf is parallel to Vi since the speed halo DM is in general much greater 
than that of the trapped DM. 


Having discussed the properties of C c and C s , we return to Eqs. (15a) and (15b). Since 


we are interested in the symmetric DM, the condition N x = N x holds. Thus, Eqs. (15a) and 


(15b) can be simplified into one, 


dN x 

dt 


= C c - C e N x + C S N X - (C a + C se )Nl 


with the solution 


where 


N x (t) = — 


C' c tanh(f/r y i) 


r a 1 - (C s - C e ) tanh(f/r A )/2’ 
1 


TA = 


(19) 


( 20 ) 


( 21 ) 


yjC c {C a + C se ) + (C 8 - c e y/A 

is the time-scale for the DM number in the Sun to reach the equilibrium, i.e., dN x /dt = 0 
The DM number reaches the equilibrium when tanh(f/rA) ~ 1. 

The DM annihilation rate in the Sun is given by 


r A = C a N x N x = C a N 2 x . (22) 

By setting C s = C se = 0, we recover the results in Refs. j44l 148] for the absence of DM 
self-interaction. By setting C e = C se = 0, we recover the result in Ref. [39], which includes 
the DM self-interaction while neglects the DM evaporation. 
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FIG. 4: The dependencies of tanh (£®/ ta ) on m x and for e 7 = 10~ 9 (left) and 10~ 10 (right). 
The regions circled by thick lines are the parameter ranges in which the DM number N XyX have 
not yet reached the equilibrium in the current epoch. The region to the left of the dash line is 
the parameter space that DM evaporation takes place. This is the non-detection region for the 
indirect DM search. The upper panel is for the isospin symmetric case while the lower one is for 
the isospin violation case. 


B. Numerical results 


To illustrate the role of the light force carrier we first identify the parameter range 
on the m x — plane in which N x (t ) has already reached to the equilibrium at the present 
epoch. In the hidden 17(1) model, isospin violation is generally introduced since r] can take 
any value [5Hj. However, to simplify our discussions, we consider two extreme scenarios in 
our calculation. One is the isospin symmetric scenario with rj — 1, and the other corresponds 
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Ey = 10 -9 


£y = 10 


-10 



FIG. 5: The theoretical predictions of cr xp with = 1 MeV (red), lOMeV (green) and 100 MeV 
(blue) are given in solid lines with e 7 = 10~ 9 (left panel) and 10~ 10 (right panel). The light and 
deep purple regions are excluded by LUX [63] for isospin symmetric and isospin violating scenarios, 
respectively. The orange band is the SIDM allowed region [221425] . See the main text for details. 


to t] = —0.7 which minimizes DM-Xenon cross section for a fixed a 


xv 


The first scenario 


corresponds to £ 7 /£z = —0.65 while the second one corresponds to £ 7 /£z = 5.65. 

In Fig. |dj we show the equilibrium regions for e 7 = 10“ 9 and 10~ 10 [60, [6T] with r\ = 1 
and r] = —0.7, respectively. The dense color regions circled by the thick lines are the 
parameter ranges where the equilibrium has not yet reached in the current epoch. We 
quantify such regions by tanh^/r^) < 0.99. To the left of the gray dashed line, DMs in 
the Sun evaporate and consequently unable to produce detectable signals. In contrast, the 
equilibrium is between the capture and the annihilation to the right of the gray dashed line. 
In this region, the captured DM number is at the maximum, hence the annihilation rate 
is at its maximum as well. It is of interest to compare the size of non-equilibrium region 
for different combinations of e 7 and in Fig. [dj One can see from Eq. (21) that a smaller 
C c leads to a larger equilibrium time scale ta provided all the other coefficients are held 
fixed. Comparing the left and right panels of Fig. [dj it is seen that the non-equilibrium 
region increases tremendously as e 7 lowering from 1CT 9 to 1CT 10 . In fact, one can see from 
Eq. © that a xp decreases by two orders of magnitude when e 7 decreases by an order 
of magnitude. Hence the DM-nucleus cross section g x a and consequently the C c are also 
significantly reduced. This leads to a larger equilibrium time scale as just argued. For 


£y = 10 , the non-equilibrium region also increases with Tj lowering from 1 to —0.7. This 
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is also due to the suppression of g x a caused by isospin violation. However, with = 1CT 10 , 
the non-equilibrium region does not change noticeably when rj changes from 1 to —0.7. For 
such an e 7 , one can show that it is C s rather than C c that dictates the equilibrium time 
scale T 4 ., i.e., the DM self-interaction dominates the capture process. 

The theoretical predictions on a xp are shown in Fig. [5] with colored solid lines. These 
predictions are according to Eq. 0 and the additional momentum transfer suppression 


given by Eq. (16). It is seen that <j xp depends on e 7 , a x and m^. However, a x is related to 
m x by a x ~ 3.3 x 10 ~ 5 (m x / GeV) from the thermal relic density of symmetric DM |3Sl Wl\. 
The regions excluded by LUX [63] for the scenarios of isospin symmetry and isospin violation 
are displayed for comparisons. The orange band is the allowed region of the SIDM parameter 
space [22D25] . 

The value of rj affects the DM-nucleus cross section as seen from Eq. 0 - Hence rj ^ 1 
(isospin violation) may weaken the direct search bound on a xp as well as suppress the capture 
rate for both the Sun and the Earth for a fixed a xp as discussed in Refs. [U3G EH] - It is clearly 
seen that the LUX bound on a xp is significantly weaken for r] = —0.7. 

The evolution behaviors of N x with different m x , and rj are shown in Fig. [ 6 j On the 
upper left panel, the difference in N x between rj — 1 and rj = —0.7 can be understood by 
the following arguments. First, the coefficient C c with rj = —0.7 is suppressed compared to 
C c with rj = 1. Hence N x with rj = 1 is greater than N x with rj = —0.7 no matter how 
large C s is. However, as m x increases, C s becomes more and more dominant over C c since 
a xx in C c grows as rn A x while <j x a in C c grows as m x (this follows from Eqs. (| 8 j) and (|9j) and 
the fact that a x grows as m x ). Hence the reduction of C c by rj = —0.7 does not produce 
significant effect on N x , given a very dominant C s . This explains why N x is less sensitive 
to rj for larger m x . Comparing upper right panel with the upper left one, we End N x is less 
sensitive to rj in the former case for the same m x . The DM-nucleus cross section a x A for the 
upper right panel is suppressed due to a smaller e 7 . Hence the effect from rj becomes less 
significant. For lower left and lower right panels of Fig. [ 6 J we take m ?J = 100 MeV. In this 
case, both g x a and a xx are suppressed. Hence the number of captured DM is much less and 
the equilibrium time scale is much larger. 
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B-y = 10“ 9 , = 30 MeV 


E y = 10~ 10 , = 30 MeV 






FIG. 6 : The N X (N X ) for e 7 = 10 ~ 9 (left column) and 1CD 10 (right column). The upper panels 
are for m^ = 30 MeV and the lower panels are for m^ = 100 MeV. Solid lines are for isospin 
symmetric case and dashed lines are for isospin violated case. Blue, green and red curves are for 
m x = 5 GeV, 20 GeV and 100 GeV, respectively. The solar age (t & ) is indicated by the gray dashed 
line. 

IV. TESTING SIDM MODEL IN ICECUBE-PINGU 

A. Neutrino signal and the atmospheric background 

Neutrino signals arise from the decays of (ft which is produced by DM annihilation, i.e., 
XX —> <ft(ft —> 4za DM annihilate to ZZ and 77 modes are suppressed by small couplings 
between y-Z and y -7 respectively. For the hidden U( 1) gauge model considered in this paper, 
<ft decays into SM particles via e 7 and mixings. Due to the kinematics constraint, e + e~ 
and neutrinos are the only decay products of (ft and the corresponding decay widths through 
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photon and Z mixings are given in Ref. [38]. The branching ratio for 0 —>■ vv is determined 
by the relative magnitudes of mixing parameters £ 1 and Ez■ Therefore it is determined by 


the parameter 7] through Eqs. (10a) and (10a). We find BR(0 —>- uv) ~ 75%, 39%, 48%, 
and 67% for 7] = 1, —0.3, —0.5, and —0.7, respectively. Furthermore, although 0 is produced 
on-shell in the solar center, it does not decay instantly but instead propagates for a certain 
distance. The lifetime of 0 is constrained by BBN, ry < 0(1) s [38], [6T]. Considering a 30 
MeV 0 with 5 GeV energy (corresponding to m x = 10 GeV), the decay time of 0 would 
be less than 170 s. One can easily estimated that, for a 0 with 170 s of decay time, it 
travels for roughly 5 x 10' km before decaying into the electron-position or the neutrino 
pair. The decay point is already outside the Sun but not yet reaching to the Earth. Given 
the distance between the Sun and the Earth at 1.5 x 10 8 km, those 0 with 30 MeV of mass 
and moving toward the Earth shall decay between Sun and Earth provided 1% is less than 
150 GeV. Hence the neutrino flux will be observed by the terrestrial detector. We have so 
far made our argument with ry = 1 s. For < 1 s, the decay point of 0 could be inside 
the Sun. In this case, the neutrino propagating distance from source to terrestrial detector 
becomes larger. However this does not affect the oscillation of neutrinos because neutrino 
propagating distances in both cases are much larger than the neutrino oscillation length. 
Hence our results on neutrino event rates are not affected. 

The argument in the last paragraph assumes 0 propagating freely inside the Sun. To 
justify this assumption, we study the interaction between 0 and hydrogen, which occurs 
through the mixing between 0 and 7 and the subsequent yp scattering. The total yp 
scattering cross section at yfs = \/T0 GeV is about 0.1 mb [68]. Taking = ICE 9 , we have 
0p total cross section as small as 10~ 46 cm 2 . With the average hydrogen number density in 
the Sun about 6 x 10 23 /cm 3 , the mean free path of 0 is much greater than the radius of the 
Sun. We note that the consideration of other chemical elements in the Sun should shorten 
the mean free path of 0 slightly. Nonetheless, this distance scale remains much greater than 
the radius of the Sun. 

Since each 0 stays at the same direction before its decay, one can view those neutrinos 
produced by 0 decays as being originated from the core of the Sun. Hence we can write the 
neutrino flux as 


d®j. 

dE „, 


= P„ 


XK) 


r A dN„ 
AttR 2 dE v 


(23) 
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with P v the neutrino oscillation during the propagation, the DM annihilation 

rate, R the distance between Sun and Earth, and dN v ./dE Vj the neutrino spectrum from 
each annihilation. The energy distribution dN v /dE v for neutrinos produced by (ft decays 
is not a simple Dirac-5 function, since (ft is highly boosted, ft has been shown that [69] 

dN 4 

-HF- = a7F 0 (^ - £-)©(£+ - Ev)> (24) 


maximum 


where 0 denotes the Heaviside step function. E± = (m x ± \J m \ ~ m|)/2 are the 
and minimum of the neutrino energy. A E = y/m^ — denotes the width of the energy 
spectrum. 

The neutrino event rate in the detector is given by 


r m x rl<f) 

N v = / —^A„(E„)dE u dn , (25) 

k dE -i 

where E t h is the detector threshold energy, A U (E U ) the detector effective area and D the 
solid angle. We study both muon track events and cascade events induced by neutrinos. The 
DeepCore detector extends the IceCube capability to probe E v down to 10 GeV mi and the 
future PINGU detector will further lower down the energy threshold E th to 0(1 ) GeV [72]. 
The angular resolution for IceCube-PINGU detector at E v = 5 GeV is roughly 10°. Hence we 
consider neutrino events arriving from the solid angle range AD = 27r(l —cos^) surrounding 
the Sun with ift = 10°. The detector effective area of IceCube is expressed as 


N a 

Ai/(E v ') Peff , r up(Ei;) V 71 i/n(E 

c.p. 


(26) 


with V e s the detector effective volume IZH [72], Na the Avogadro number, M ice the molar 
mass of ice, n Pi „ the number density of proton/neutron per mole of ice, and <J up , n the neutrino- 
proton/neutron cross section. One simply makes the replacement v —> v for anti-neutrino. 


The atmospheric background is similar to Eq. (25), by replacing d§ v fdE v with atmo¬ 
spheric neutrino flux, 


A^atm 


d® 


atm 


dE v 


-A u (E u )dE u dQ. 


The dAO m /dE u is taken from Ref. 
signal. 


We set E max = m x in order to compare with DM 
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Track 


Cascade 




FIG. 7: The IceCube-PINGU sensitivities to for different values of r\. The left panel is for 
track events and the right panel is for cascade events. The sensitivity curves and LUX excluded 
regions correspond to e 7 = 10~ 9 . The results corresponding to e 7 = 1CU 10 are not shown since 
BBN constraint rules out most of the parameter space plotted here. 

Track Cascade 


£max [GeV] N* tm N »M IV* tm iV° M 


5 

7146 

76 

9874 

90 

10 

10280 

91 

13775 

105 

50 

21680 

132 

21803 

132 

70 

23584 

138 

23111 

136 

100 

26610 

146 

24363 

140 


TABLE I: The annual signal and background event numbers for reaching 2cr detection significance 
in 5 years. 

B. The IceCube sensitivity to the mass of the force carrier 


We present the sensitivity as a 2cr detection significance in 5 years, calculated with the 
formula, 


s 

Vs + b 


2 . 0 , 


where s is the DM signal, b is the atmospheric background, and 2.0 is referring to the 
2cr detection significance. The detector threshold energy E t h is set to be 1 GeV, which 
will be achieved by the future PINGU detector. However, for E u > 10 GeV, we adopt the 
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DeepCore effective area in our calculation. The annual signal and background event numbers 
for reaching 2 a detection significance in 5 years are listed in Table |T| We present the track 
and cascade event numbers separately. 

In Fig. [7j we present the IceCube-PINGU sensitivities to with different values of 77 in 
colored solid lines. Regions above these lines are beyond the reach of the detector within 5 
years for a 2 cr detection significance. Furthermore, each sensitivity curve terminates at the 
evaporation mass scale below which the DM signature from the Sun is suppressed. The left 
panel is for track events and the right one is for cascade events. BBN excludes < 20 MeV 
for e 7 = 10 _1 °, while it excludes < 0.3 MeV for e 7 = 10~ 9 . For simplicity, we only present 
results for e 7 = 10 -9 since BBN constraint for e 7 = 10 -1 ° rules out most of the m^ parameter 
space shown here. The orange band is SIDM allowed region. Those color shaded regions are 
excluded by LUX for different 77 values. The gray dashed lines indicate the relation between 
a x and m x required by the thermal relic density. 

It is of interest to compare results with different 77 values. It is seen that the sensitivity 
curve with 77 = 1 differs significantly from those with other 77 values. However, sensitivity 
curves with —0.3 < 77 < —0.7 do not differ much from one another. Clearly C s rather than 
C c dominates the DM capture in this range of 77 . 

Comparing left and right panels of Fig. [7J one can see that cascade events can probe a 
larger part of m x — parameter space. Furthermore, only cascade events can probe into 
the SIDM allowed region. Comparing IceCube-PINGU sensitivity with LUX constraint, the 
former complements the latter for low-mass DM until the evaporation mass scale m x < 
4 GeV. Such a complementarity also occurs due to isospin violation. It is well known 
that the direct search bound on cr xp could be significantly weaken by isospin violation. In 
fact, one can see that the LUX excluded region shrinks as 77 decreases from 1 to —0.7. On 
the other hand, the indirect DM signature from the Sun is less affected by 77 . Hence the 
complementarity of two searches becomes more apparent when a xp is subject to more severe 
isospin violation suppression. 

V. SUMMARY 

We have proposed to test SIDM model with hidden 1/(1) gauge symmetry by searching 
for DM-induced neutrino signature from the Sun. The SIDM model under consideration 
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gives rise to DM-nucleus scattering, DM-DM elastic scatterings, and DM-DM annihilation. 
These three processes determine the neutrino flux resulting from DM annihilation in the 
Sun. We have presented the IceCube-PINGU sensitivities to the parameter space of SIDM 
model. We compare these sensitivities to existing constraints set by LUX experiment. We 
have shown that the indirect search complements the direct one in two ways. First, the 
indirect search is more sensitive to light DM in GeV mass range. Second, the direct search 
constraint on SIDM parameter space can be significantly weaken by isospin violation while 
the indirect search sensitivity is less affected by this effect. 
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